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Abstract 

We prove that the binary complexity of solving ordinary polynomial differential equa- 
tions in terms of Puiseux series is single exponential in the number of terms in the series. 
Such a bound was given by Grigoriev [TO] for Riccatti differential polynomials associated 
to ordinary linear differential operators. In this paper, we get the same bound for ar- 
bitrary differential polynomials. The algorithm is based on a differential version of the 
Newton-Puiseux procedure for algebraic equations. 

Introduction 

In this paper we are interested in solving ordinary polynomial differential equations. For 
such equations we are looking to compute solutions in the set of formal power series and 
more generality in the set of Puiseux series. There is no algorithm which decide whether 
a polynomial differential equation has Puiseux series as solutions. We get an algorithm 
which computes a finite extension of the ground field which generates the coefficients of the 
solutions. Algorithms which estimate the coefficients of the solutions are given in |12} I13j. 

In order to analyse the binary complexity of factoring ordinary linear differential op- 
erators, Grigoriev [10] describes an algorithm which computes a fundamental system of 
solutions of the Riccatti equation associated to an ordinary linear differential operator. 
The binary complexity of this algorithm is single exponential in the order n of the linear 
differential operator. There are also algorithms for computing series solutions with real 
exponents [HI [IJ [2] and complex exponents [8]. 

Let K = Q(Ti, . . . ,Ti)[rj\ be a finite extension of a finitely generated field over Q. 
The variables T\ , . . . , Tj are algebraically independent over Q and r] is an algebraic element 
over the field Q(Ti, • • • , Tj) with the minimal polynomial <p G Z[Ti, . . . ,Tj] [Z]. Let K be an 
algebraic closure of K and consider the two fields: 

L = \J v&i .K((x$)), C = U uen *K((x^)) 

which are the fields of fraction-power series of x over K (respectively K), i.e., the fields of 
Puiseux series of x with coefficients in K (respectively K). Each element ip S L (respectively 
ip € C) can be represented in the form tp = ^2 ie qCiX l , C{ € K (respectively q € K). The 
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order of ip is defined by ord(tp) := min{i £ Q,Cj / 0}. The fields L and £ are differential 
fields with the differential operator 

ieQ 

Let F(uq, . . . , y n ) be a polynomial on the variables yo,---,y n with coefficients in L and consider 
the associated ordinary differential equation F(y, , • • • , ^r) = which will be denoted by 
= 0. We will describe all the solutions of the differential equation F(y) = in C by a 
differential version of the Newton polygon process. First write F in the form: 

F= £ haX^-'-y^, fi, a €K 

where a = (ao, • • • , a n ) belongs to a finite subset A of N n+1 . The order of F is defined by 
ord(F) := min{i € Q; f^ a / for a certain a}. Without loss of generality we can suppose 
that each coefficient /j iQ 6 Z[Ti, . . . , Tj][r/] and so it can be written in the form 

fi,a= ^2 ^j,h,-,3l^l' ' ' ' Ff l rf , bjj lj ,,,j l £ Z. 
3,31,— ,31 

We define the degree of F w.r.t. x by deg^F) = max{i G Q; fi, a ^ for a certain a} 
(it can be equal to +oo), the degree of F w.r.t. Ti,...,Ti by deg Tl T (F) = 
max{deg Tlj iTi (/j jQ ,); i € Q, a € A}. We denote by 1(b) the binary length of an inte- 
ger b. The binary length of F is defined by 1(F) = max{I(/j if( ); i £ Q, a £ ^4} where l(fi t a) is 
the maximum of the binary lengths of its coefficients in Z. We can define in the same manner 
the degrees and the binary length of <f>. To estimate the binary complexity of the algorithm 
of this paper we suppose that deg^(^) < do, deg Tl T t {4>) ^ ^ deg yo ^ yn (F) < 

d, deg Tli < c?2, deg 2 ,(i ? ) < c?3 (cfo can be equal to +00) and 1(F) < M 2 . 

In the sequel, we will use the following quantities: for each ieQ, let 

R(i) = d 2 (dd d l )°^ 

and 

S(i) = n i d 2 (dd d 1 ) o(a \M 1 M 2 )°^ log 2 (dd 3 ). 
We can now describe the main theorem of this paper: 

Theorem 0.1 Let F(y) = be a polynomial differential equation with the above bounds. 
There is an algorithm which computes all Puiseux series solutions of F(y) = with coeffi- 
cients in K , i.e., all solutions of F(y) = in C. Namely, for each solution ip = J2ieQ c « x * ^ ^ 
of F(y) = 0, the algorithm computes an integer v E N* such that for each i G Q, it computes 
a finite extension K\ = K[8] of K where 9 is an algebraic element over K computed with 
its minimal polynomial $ £ K[Z] such that J2j<i j^Q c j x ^ £ FJ\((x~)). Moreover, for any 
j < i, j £ Q, we have the following bounds: 

- deg z ($) < d\ 
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- deg Tl) ... iTi ($), deg Tli ... )T( (cj) < R(i). 
-l($),l( Cj )<S(i). 

- The binary complexity of this computation is S(i). 

Remark 0.2 i) By the corollary of Lemma 3.1 of 11 Jj /. the integer v of Theorem \ 0.1\ is 
constant, i.e., independent of i. This constant depends only on the solution tjj. 

ii) In general, we cannot compute a finite extension K\ of K which contains all the 
coefficients of all the solutions of F(y) = in C either an integer v G N* such that all the 
solutions of F{y) = (in C) are in Ki((x*)). Namely, if we consider the polynomial 

F{yo, 2/1,2/2) = xy y 2 - xy\ + y yi 
then if) = cx^ is a solution of F(y) = in C for all c G C and all fj, G Q. 

In a forthcoming paper we try to get polynomial binary complexity for computing Puiseux 
series solutions of polynomial differential equations as in the algebraic case [5j [T] . 

The paper is organized as follow. In section 1, we give a description of the Newton 
polygon associated to polynomial differential equations. Operations on these Newton 
polygons will be discussed in section 2. The algorithm with its binary complexity analysis 
are described in section 3. 

1 Newton polygons of polynomial differential equations 

Let F be a differential polynomial as in the introduction. We now define the Newton polygon 
of F. For every couple (i,a) G Q x A such that fi <a 7^ (i.e., every existing term in F) we 
mark the point 

Pi,a '■= {i — ai — 2«2 - • • • - na n , a + a\ H h a n ) G Q x N. 

We denote by P(F) the set of all the points Pi )Q ,. The convex hull of these points and (+00, 0) 
in the plane M. 2 is denoted by M(F) and is called the Newton polygon of the differential 
equation F(y) = in the neighborhood of x = 0. If deg yo ^ yn (F) = m, then M(F) is 
situated between the two lines y = and y = m. For each (o, b) G Q 2 \ {(0, 0)} we define the 
set 

N(F, a, b) := {{u, v) G P(F), V(«', v') G P(F), av! + bv' >au + bv}. 

A point P it a G P(F) is a vertex of the Newton polygon Af(F) if there exist (a, b) G Q 2 \{(0, 0)} 
such that N(F,a,b) = {Pi y0l }. We remark that J\f(F) has a finite number of vertices. A pair 
of different vertices e = (Pj, a , Pi\ a ') forms an edge of N(F) if there exist (a, b) G Q 2 \ {(0, 0)} 
such that e C N(F,a,b). We denote by E(F) (respectively V(F)) the set of all the edges e 
(respectively all the vertices p) of Af(F) for which a > and b > in the previous definitions. 
It is easy to prove that if e G E(F), then there exists a unique pair (a(e),6(e)) G 1? such 
that GCD(a(e),b(e)) = 1, a(e) > 0, 6(e) > and e C N(F, a(e),b(e)). By the inclination of 
a line we mean the negative inverse of its geometric slope. If e G E(F), we can prove that 
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the fraction p e := G Q is the inclination of the straight line passing through the edge 

e. If p G V(F) and N(F,a,b) = {p} for a certain (a, b), then the fraction p := ~ G Q is the 
inclination of a straight line which intersects M{F) exactly in the vertex p. 

For each e € E{F) we define the univariate polynomial (in a new variable C) 

ff(F,e)(C) := £ / iiQ C a »+ ai +-+ a «( Me )r • • • GuX 1 G # [C], 

P< l0 eiV(^,o(e),6(e)) 

where (p e )k '■= pe(pe — 1) • • • (p e — k + 1) for any positive integer fc. We call H^ Fe ^{C) 
the characteristic polynomial of F associated to the edge e G E{F). Its degree is at most 

"» = de Sy ,...,y n ( F ) ^ d - 

If if) G £ is a solution of the differential equation F{y) = such that ord(ip) = [i e , 
i.e., V has the form ip = ^ieQ i>fJie Qx 4 , Cj G if, then we have that H^ F ^{c^ e ) = 0, i.e. c Me is 
a root of the polynomial Hr Fe \ in K. This condition is called a necessary initial condition to 
have a solution of F(y) = in the form of ?/> (see Lemma 1 of [2]). In fact, Hip^ic^) equals 
the coefficient of the lowest term in the expansion of F(ip(x)) with indeterminates p e and 
c Me . Let A (Fje) ■= { C £ K,c^ 0,H {F ^ e) (c) = 0}. 

For each p = (u,v) G V(F), let p\ < p2 be the inclinations of the adjacent edges at 
p in Af(F), it is easy to prove that for all rational number p = a G N*, 6 G N such that 
N(F, a,b) = {p}, we have p\ < p < p2- We associate to p the polynomial 

which is called the indicial polynomial of F associated to the vertex p (here p is considered 
as an indeterminate). Let Hm p \(C) = C v h^ F ^(p) defined as above for edges e G E{F). Let 

a (f, p ) ■= e Q, mi < /i < h {F , p) (p) = 0}. 

Remark 1.1 Xei p = (u, i>) G ^(i 7 ) anc? e be the edge of M{F) descending from p, then 
^(F,p)(^e) is the coefficient of the monomial C v in the expansion of the characteristic polyno- 
mial of F associated to e. 

2 Some operations on Newton polygons of differential poly- 
nomials 

2.1 Relation between Newton polygons of a differential polynomial and its 
partial derivatives 

Let F be a differential polynomial as in the introduction. Write F in the form 

F = F + ■ ■ ■ + F d 

where F s = ^igq \a\=s /i,a^?/o ' ' ' Vn n is the homogeneous part of F of degree s with respect 
to the indeterminates yo, . . . , y n , a = (otQ, . . . , a n ) G A and \a\ = ao + • • • + a n is the norm of 
a. Then the ordinate of any point of P(F S ) is equal to s and 

P(F) = U < s < d P(F s ). 
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Let < j < n. If there exists an integer k > 1 such that for all 1 < s < d (such that F s ^ 0), 
D s j := deg r (F s ) > k then we can prove that P is a translation of -P(F) defined by the 

point (kj, —k), i.e., 

'd^Fx r fd k F^ 



For any (a, 6) G Q 2 \ {(0, 0)}, we have 

^(0, a, &) = N{F, a, b) + {(kj, -k)}. 

Thus the edges of TV (fjpr) are exactly the translation of those of Af(F). For each e G E^^^j, 
its characteristic polynomial is 

H (9kF AC) = £ ^ aC *o + a 1+ -+a B - ^a,-)*^)? 1 • • • foe)?"* ■ ■ ■ (rf € 

9 4 ' P i}a eN{F,a(e),b(e)) 

For each p G ^(^f")' * ts ^dicia/ polynomial is 



/l (e^F p )(^)= X] fiA a 3)k{Ve)T ■ ■ ■ (Me)? • • • (Ai e )£" € -K"[/i]. 

Let < ji 7^ j2 < n. If there exist integers fci, > 1 such that for all 1 < s < d, D s j 2 > k 2 

aIld deg %! (^f) ^ ^ then 

fcl fc2 ) = W + {(hji + — fci - k 2 )} 
Uj 1 Uj 2 

and then 

/ f)ki+k 2 

M {-]fTrT7) = W F ) + {(^'i + ^J2, -fci - k 2 )}. 
Vj 1 Uj 2 

For any (a, 6) G Q 2 \ {(0, 0)}, we have 

,Qki+k 2 

N( ki k2 ) = N(F, a, b) + {(hji + k 2 j 2 , -h - k 2 )}. 
For each e G E 1 ( 9 ^ 2 fc f ) , its characteristic polynomial is H, g k 1 +k 2F x(C) = 



£ /^C' a0+Q1+ -" +Q "" fel " fe2 K 1 )fe 1 (^)fe 2 (^)r • • • (Me)? 1 ^ 1 • • • (Me)?^ 2 • • • 

Pi,«eAT(F,a(e),&(e)) 

Lemma 2.1 Zei k > 1 be an integer such that for all 1 < s < d (such that F s ^ 0) and for 
all < j < n we have D s j > k. For any e G E(F), the k-th derivative of H^ Fe ^ G K[C] is 
given by the formula 

H (t)^ = £ Mi ■ ■ ■ i^t n H { dkF e) (C). 

0<fc ,...,fc„<n dyo°--Syt-' 

where the sum ranges over all the partitions (/cq, • • • , k n ) of k, i.e., ko + • ■ ■ + k n = k. 
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Proof. By induction on k. □ 



2.2 Newton polygons of sums of differential polynomials 

Let F and G be two differential polynomials of degree less or equal than d. Write F = 
Fq + • • • + Fd and G = Gq + • • • + Gd where F s (respectively G s ) is the homogeneous part 
of F (respectively G) of degree s with respect to the indeterminates yo, ■ ■ ■ ,y n - F° r each 
< s < d such that F s + G s / 0, let P s = (s\,s) be the point of the plane defined by 
s\ := min{u; (u,s) G P(F S + G s )}. Then it is easy to prove that the convex hull of all 
the points P s for all < s < d and (+oo, 0) in the plane M 2 is the Newton polygon of the 
differential polynomial F + G. 

2.3 Newton polygons of evaluations of differential polynomials 

Let F be a differential polynomial as in the introduction, ^ c £ K, /i G Q and 
G(y) = F(cx M + y). We will discuss the construction of the Newton polygon of the differential 
equation G(y) = for different values of c and \x. 

For each differential monomial m(y) = f\,aX l y^ ■ ■ ■ y%™ of F with corresponding point 
p, compute 



Remark that the corresponding points of the differential monomials of m(cx^ + y) have ordi- 
nate less or equal than s and lie in the line passing through p with inclination fx. There are 
two possibilities for \i: 

Lemma 2.2 If ji = ji e is the inclination of an edge e G E(F). For any < s < d, the vertex 
ofJ\f(G) of ordinate s corresponds to the differential monomial of G with coefficient equals to 



where the sum ranges over all the partitions (ko, . . . , k n ) of s. Its x-coordinate is the minimum 
of the quantities i + /u(ao + ■ ■ ■ + a n — s) — a± — 2«2 — • • • — na n for i G Q and a G A. 
If Hi < fj, < fj,2 where fj,i and [12 are the inclinations of the two adjacent edges of a vertex 
p = (u,v) G V(F), then for any < s < d, the vertex of N{G) of ordinate s corresponds to 
the differential monomial of G with coefficient equals to 



m{cx» + y) = fi, a x l {cx» + y ) ao (c^- 



+ yi )^...(c(ti) n x»- n + y n ) 




0<fco<- • -<k„<n 





C' 



,v—s 



0<ko<---<k n <n 





where the sum ranges over all the partitions (ko, . . . , k n ) of s. 



Proof. Let h = \x e for e G E(F) and compute 



G(y)= Yl fi,aX i (cx^ + y r o (c^- 1 +y 1 ) a ^--(c^) 



+ y n 



) 
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For each < s < d, compute G s the homogeneous part of G of degree s in yo,...,y n . We 
remark that for fixed i and a, all the differential monomials of G s have the same corresponding 
point which is 

(i + ju(«o + ■ ■ ■ + a n — s) — a\ — 2a 2 — • • • — na n , s). 

The x-coordinate of the vertex of M(G) of ordinate s is the minimum of the x-coordinates 
i + /i(«o + ■ ■ ■ + a n — s) — a\ — 2a 2 — • • • — na n for i G Q and a £ A. This minimum is 
realized by the points Pj )Q , G iV(-F, a(e), 6(e)). This proves the lemma taking into account the 
formula for the characteristic polynomial of the derivatives of F in the previous subsections. □ 

The following lemma is a generalization of Lemma 2.2 of [10J which deals with the 
Newton polygon of the Riccatti equation associated to a linear ordinary differential equation. 

Corollary 2.3 Let \i = fi e be the inclination of an edge e G E(F). The edges of J\f(G), 
situated above the edge e are the same as in M(F). Let an integer si > be such that 
q s (c, fj> e ) = for all < s < s\ and q si (c, fi e ) ^ then Af(G) has a vertex of ordinate si and 
it has at least s\ edges with inclination greater than \i e . 

Proof. Let p Sl be the vertex of M(G) of ordinate si and x-coordinate x Psi the minimum of 
the values i + fi(ao + ■ ■ ■ + a n — s\) — a\ — 2«2 — • • • — na n for i S Q and a € A (by Lemma l2T2j) . 
We have q Sl -i(c, fi e ) = 0, then the x-coordinate of the vertex p Sl -i of ordinate s\ — 1 is strictly 
less than the minimum of the values i + fi(ao + • • • + a n — s\ + 1) — ot\ — 2a 2 — • • • — na n for 
i G Q and a G A. Thus the inclination of the edge joining p Sl and p s ±-i is greater than [i. □ 

Corollary 2.4 Let /x = [i e be the inclination of an edge e G E{F). If Hm e \(c) = then the 
intersection point of the straight line passing through e with the x-axis is not a vertex of N{G) 
andN{G) has an edge with inclination greater than /i e . 

Proof. We have qo(c, /j, e ) = H^ F ^{c) = 0, then s\ > 1. This proves the corollary by applying 
Corollary [231 □ 

3 Differential version of the Newton-Puiseux algorithm 

We describe now a differential version of the Newton-Puiseux algorithm to give formal 
Puiseux series solutions of the differential equation F(y) = 0. The input of the algorithm 
is a differential polynomial equation F(y) = with the bounds described in the introduc- 
tion. The algorithm will construct a tree T which depends only on F and on the field 
K. The root of T is denoted by tq. For each node r of T, it constructs the following elements: 

- The field K T which is a finite extension of K. 

- The primitive element 6 T of the extension K T of K with its minimal polynomial (fi T G 
K[Z]. 

- An element c T G K T , a number fi T G Q U {— oo, +00} and an element y T = c T x^ T + y T1 G 

1 

Kt^x 17 ^)) where r is a descendant of n (here /j, t > ji Tl ) and u(t) G N*. 

1 

- The differential polynomial F T {y) = F{y + y T ) with coefficients in K T ((x v ^)). 

We define the degree of r by deg(r) = /i T e Q, we have deg(r) = deg x (y T ) if c T 7^ 0. 



7 



The level of the node r, denoted by lev(r), is the distance from to to r. 

For the root To we have K T0 = K, 8 TQ = 1, (p TQ = Z — 1, c T0 = y T0 = 0, deg(ro) = 
/"to = -oo,z/(t ) = 1 and F TQ (y) = F(y). 

A node r of the tree T is a leaf of T if for each e € E(F T ) and for each p £ ^(-^V) 
we have fi e < deg(r) and fi2 < deg(r) and y = is a solution of -F T (y) = 0, where < /12 
are the inclinations of the adjacent edges at p in Af(F). 

The algorithm constructs the tree T by induction on the level of its nodes. We sup- 
pose by induction on i that all the nodes of T of level < i are constructed. Denote by % the 
set of these nodes. At the (i + l)-th step of the induction, for each node r of level i which is 
not a leaf of T we consider the following two sets: 

- E'{F T ) = {e € E(F T ), He > deg(r)} and 

- V'{F T ) = {pe V(F T ), V2 > deg(r)}. 

For each e E E'(F T ), compute a factorization of the polynomial H(p T ^(C) £ K T [C] 
into irreducible factors over the field K T = K[9 T ] in the form 

H {FT ^{C) = X e \{Hf 
j 

where ^ \ e G K T , kj € N* and Hj € -ftT T [C] are monic and irreducibles over K T . We can do 
this factorization by the algorithm of [31111 El [9]. The elements of the set A^p T ^ correspond 
to the roots of the factors Hj ^ C. We consider a root Cj £ K for each factor Hj ^ C and 
we compute a primitive element #j, e ,T of the finite extension K T [cj] = K[9 T , Cj] of K with its 
minimal polynomial 4>j, e ,r € K[Z] using the algorithm of [HJ El [9] . 

For each root Cj of Hj ^ C we correspond a son a of r such that CT = #j, e ,T> the 
field = /f[#j ieiT ] and the minimal polynomial of B a over K is (f) a = (f>j, e ,r- More- 
over, c a = Cj, = n e , y a = c a x^ + y T and F a (y) = F(y + y a ). For i/(cr), we take 
v(cr) = LCM(u(T),a(e)) for example. 

For each p € V'(F T ), we consider the indicial polynomial h(^p Tp ^{p,) € if T [/i] of F T as- 
sociated to p. To each € A^ FtP ^ such that /x > deg(r) and 0/c£K (where c is given by 
its minimal polynomial over K), we correspond a son a of r such that 9 a = c, c CT = c, z/ CT = /x. 
This completes the description of all the sons of the node r of the tree 7~. 

Remark 3.1 i) If (E'(F T ) / or V'(F T ) ^ ®) and y = is a solution of F T (y) = then one 
of the sons of r is a leaf a for which F a = F T , \i a = +oo and c a = 0. 

ii) For any node t of T such that deg(r) ^ oo, if y = is not a solution of F r (y) = then 
E'{F T ) / by Corollary 

3.1 Determination of the solutions of = in £ by the leaves of T 

Let £Y be the set of all the vertices r of T such that either deg(r) = +oo and for the ancestor 
t\ of r it holds deg(ri) < +oo or deg(r) < +oo and r is a leaf of T. For each t £ U, there 



8 



exists a sequence (rj(r))j>o of vertices of T such that tq(t) = tq and Tj + i(r) is a son of u(t) 
for all i > 0. For each t€W, the element 

yr = J^c^a^iM G K T ((x^)) 

is a solution of -F(y) = 0. In fact, there are two possibilities to r: if deg(r) = +00 then y = 
is a solution of F T1 (y) = where r is a son of t\ and so y T1 is a solution of = 0. If 

deg(r) < +00 and r is a leaf of T then y = is a solution of F T (y) = F(y + y r ) = and so 
y T is a solution of F{y) = 0. This defines a bijection between U and the set of the solutions 
of F(y) = in C. 

3.2 Binary complexity analysis of the Newton-Puiseux algorithm 

We begin by estimating the binary complexity of computing all the sons of the root To of 
T. For each e G E(F), we consider the polynomial Htp^ e \{C) G its degree w.r.t. C 

(respectively T\, . . . ,3]) is bounded by d (respectively d 2 ). We have fi e < and its binary 
length is Z(/U e ) < 0(log 2 (dd3)) (using the fact that \x e is the inclination of the straight line 
passing through e). Then the binary length of Hi Fe \{C) is bounded by M 2 + ndO{\og 2 {dd^)). 
By the algorithm of [HJ HJ El [9] the binary complexity of factoring H( F ^(C) into irreducible 
polynomials over K is 

(dd 1 d 2 )°«(n ( i M 1 M 2 log 2 (d ( i3)) O(1) . 

Moreover, each factor Hj G K[C] of Htp je \(C) satisfies the following bounds (see Lemma 1.3 
of 0): deg c (i2j) < d, deg Tl ^ m {Hj) < d 2 (cta di)°W and 

Z(i3j) < n/d 2 (dd di) o(1) MiM 2 log 2 (dd 3 ). 

By the induction we suppose that the following bounds hold at the i-th step of the algorithm 
for each node r of T of level i: 

- deg z (<f> T ) < d\ 

- deg Tli ... jTi ((/> T ), deg Tl . jT; (c T ) < R(i). 

- 1(4> T ), l(c T ) < S(i) where R(i) and S(i) are as in the introduction. 

- fJ-r < Kit) an( i tnen K^t) < 0(log 2 (idd 3 )). 

Then we have the following bounds for the differential polynomial F T {y) = F{y + y T ) G 
K T {(x^)))[y ,...,y n ]: 

- deg yo _ yn (F T ) < d. 

- deg Tl) ... )Tl (F T ) < d 2 + ddeg Tlj ... )Tj (c r ) < R(i). 

- deg x (F T ) <d 3 + dn T <(i + l)d 3 . 
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- l(F T ) < M 2 + dl(c T ) < S(i). 

We compute a primitive element rji of the finite extension K T over the field Q(Ji, . . . , T;), 
i.e., K T = K[9 T ] = ®(T 1 ,...,Ti)[rj\[0 T ] = Q(Ti, . . . ,T,)foi] by the corollary of Proposi- 
tion 1.4 of [10] (see also section 3 of chapter 1 of [6j). Moreover, rji = rj + 7# T where 
< 7 < [K T : Q(Ti, . . . , T;)] = deg z (4>) deg z (cj) T ) < d l do and we can compute the monic 
minimal polynomial (j>\ € Q(Ti, . . . , of 771 which satisfies the following bounds: 

- deg z Oi) < cP d 
-deg Tli ... iTi (<Ai)<(dMo)° (1) 

- K<f>i) < S(i). 

- This computation can be done with binary complexity S(i). 

For each e G E'(F T ), we consider the polynomial H( Fr ^(C) € K T [C], its degree w.r.t. 
C (respectively Ti, . . . ,Ti) is bounded by d (respectively R(i))- We have /z e < (i + l)(-#) 
and its binary length is l(/J, e ) < 0^1og 2 ((i + 1)^3)J- Then the binary length of H^ Ft ^(C) 
is bounded by 

l(F T ) + ndl(fi e ) < S(i). 

By the algorithm of [3j SI El [9] the binary complexity of factoring Hrp T>e \(C) into irreducible 
polynomials over K T = Q(T%, . . . , Ti)[rji] is S(i). Moreover, each factor Hj G K T [C] of 
H( Ft ^(C) satisfies the following bounds: 

- degdHj) < d 

- degT^iHj) < R(i). 

- l(Hj) < S(i). 

Let Cj G K be a root of Hj. We can compute by the corollary of Proposition 1.4 
of [10] a primitive element %,e,r of the finite extension -fCj-fcj] = K[9 T ,Cj] of K with its 
minimal polynomial (j)j,e,T € K[Z]. We can express 0^ = Qj e ,T i n the form = r + TjCj 
where < 7^ < deg z (<p T ) deg c (Hj) < d t+l and c a = Cj in the form 

Ca = b t et 

0<<<d 4 + 1 

where bt G K. Moreover, the following bounds hold: 

- deg z (^) < d i+l . 

- degTi,... ,1,0a), deg Tli ... >Tl (b t ) < R(i). 

- KK), Kh) < S(i). 
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-/V = Me< + and then/(/v) < o(log 2 ((z + l)cM 3 )) . 

- This computation can be done with binary complexity S(i) and thus the total bi- 
nary complexity of computing all the sons a of r is S(i). 
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